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ABSTRACT 

The nonlinear space-time neutron flux equation with 
negative prompt feedback and one-group delayed neutron is 
reduced by the use of a nonlinear transformation to a par- 
tial differential equation, in which the nonlinear term 
represents a small perturbation. The general procedure of 
solution for the resulting weakly nonlinear initial-boundary- 
value problem is then established by means of the method of 
successive approximation. Convergence of the analytical 
solution is also discussed. The solutions to a slab reactor 
core and a cylindrical reactor core are investigated here. 
Asymptotic stable equilibrium states are derived from each 
of these solutions. The present results are consistent with 

those obtained from previous stability analysis for the 

2 

generalized buckling K greater or less than y Q ,.the lowest 
eigenvalue of the associated linear Helmholtz equation. 
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NOMENCLATURE 



r 

x 

t 

V 

D 

k 

k„ 

k o 

K 

C 

C P 

T 

T c 

T 

o 

e 

h 

C o 

D 

A 

V- 

a 

o 

3 

E 

a 

E f 



position vector 
dimensionless position vector 
time 

neutron speed 
diffusion coefficient 
multiplication factor 

multiplication factor of an infinite reactor 
multiplication factor at steady state 
generalized buckling 

concentration of one-group precursor 

specific heat of fuel 

temperature 

temperature of coolant 
temperature at steady state 
average energy released per fission 
convection heat transfer coefficient 

concentration of one-group precursor at steady state 

domain of a reactor core 

surface area of a fuel element 

volume of a fuel element 

feedback coefficient 

fraction of one-group delayed neutron 
macroscopic absorption cross section 
macroscopic fission cross section 
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T dimensionless time 

<f> neutron flux distribution 

<f> neutron flux at steady state 

X decay constant of one-group precursor 

p fuel density 

^ neutron flux difference, <f> - <f> 

£ precursor concentration difference, C - C q 

0 temperature difference, T - T q 

Vi, v Eigenvalues 
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I. 



INTRODUCTION 



The stability of nonlinear reactor kinetics has been 
investigated by several authors. Kastenberg and Chambre [1] 
have established stability theorems for the nonlinear space- 
time reactor kinetics equations by the method of comparison 
functions. Garabedian and Lynch [2] and Scalettar [3] have 
treated the nonlinear equations by modal expansions in 
space and then solved the resulting ordinary differential 
equations numerically. Kaplan, Margolis and Harris [4] have 
applied a pertubation technique as well as a modal expansion 
in space to obtain approximate solutions for nonseparable 
transients. Other numerical techniques were summarized by 
Kaplan, Henry, Margolis and Taylor [5]. 

In this thesis , the reduction of the nonlinear initial- 
boundary-value problem to a weakly nonlinear one is estab- 
lished by a nonlinear transformation. The transformed 
equation is now solved analytically by means of the method 
of successive approximations. As an application of this 
method, the asymptotic stable equilibrium states are 
discussed for an infinite slab reactor core and a finite 
cylindrical core. 
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II. DERIVATION OF THE SYSTEM EQUATIONS 



A. GENERAL EQUATIONS 

In this thesis, the reactor is assumed to be an one- 
velocity, bare, homogeneous and Newtonian cooling system 
with one group of delayed neutrons. The equations describ- 
ing the neutron and energy distributions and the concen- 
tration for one-group delayed neutron precursor are given 
by Meghreblian and Holmes [6] : 

= DV Z <P(C,-t) -h £ a [(i-p)k-']0(^t) + XC(C.-t) (2.1) 

_ _\C(r,-t) + fit, $(£'*) ( 2 . 2 ) 

Bt' 

(cfek.-t.) = [T(r,t)-T c (c,<-)] (2.3) 

where the resonance-escape and fast nonleakage probabilities 
to thermal are assumed to be unity. The boundary conditions 
of equations (2.1), (2.2) and (2.3) are given by: 

= C{C Sl t) =T(r 51 -t) = o -for -t> o (2.4) 

and the initial conditions are 

0(r,o) = -f(r) ; C(r , 0 ) - q{£) ■ T (r ,o) = 4uc) (2.5) 

where indicates all points on the boundary of the reactor. 

The nonlinear effects on the kinetics of a nuclear 
reactor are usually described by a single parameter the 

feedback coefficient of the reactor. The effects of 
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temperature on the multiplication factor are more pronounced 
than the effects of temperature on the other reactor 
parameters. Small changes in temperature cause small changes 
in all the reactor parameters. However, the multiplication 
factor is very close to unity so that the quantity, 

(l-g)k-l is very sensitive to even small changes in k. Hence, 
it is logical to assume that all the nonlinear coupling 
appears through this multiplication factor and the feedback 
equation is given by 



The dependence of the multiplication factor on the 
temperature and the appearance of neutron flux <f> in equations 
(2.2) and (2.3) make the system defined by equations (2.1), 
(2.2) and (2.3) intrinsically nonlinear and coupled. 
Introducing the new dimensionless variables 



. 4(T) = k 0 + [T(r.-0-T o (r ) ] 



( 2 . 6 ) 



T = v£ a t 



(2.7) 




the above system of equations becomes: 




( 2 . 8 ) 



ac(x, t) = _ _2lco,~'> -i_ £A<K*' T ) 



(2.9) 




ilT(i.t) =. 




( 2 . 10 ) 
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Defining the constants: 



y — dJh . y _ 6 If • y _ . 

°\ v ; c z - pc p \/£ r ; s 3 — pco^V-T- ' A c — 









( 2 . 11 ) 



equations (2.8), (2.9) and (2.10) can be written as follows: 

'!!(£.'?) - + + wl c C(*.t) (2.12) 









(2.13) 



= w L*,**) - n [tcx,t) -r c a)] 



(2.14) 



9 cb 9 C 3T 

From these equations , putting = 0 , - 5 — = 0 and * — = 0 

d T d T d T 

yields the steady state equations: 



o - ? 2 <&0O + [(»-f)fc.-«]£(*) + v A c C a ci) 



(2.15) 



0 =-\C e C ^ + V, 



(2.16) 



0 - - * 3 [t o (£)-T c C*)] 



(2.17) 



which are assumed to be true for x < 0 . 

Subtracting equation (2.15), (2.16) and (2.17) from 

equations (2.12), (2.13) and (2.14) respectively yields: 






(2.18) 



V L 



- -f. r, 4>cz,'? ) 



(2.19) 






( 2 . 20 ) 
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where 



0-p)K- 1 - o-pj-fe - > 



and ’ a (z) 

>(*,T> = CM-C # ^) (2.21) 

■&(JL,T) = Tu.t) - T 0 (£) 

and k(T) is defined by the feedback equation (2.6) 

•&(t) = 4 (2 * 22) 

The boundary conditions are given by : 

<y(£,/t)= M-j/'O = •0’(£*,'r) =° (2.23) 

where x designates the position vector on the boundary of 
s 

the reactor and the initial conditions are: 

y (*,<>) = R*) 

$(*,<>) = o (2.24) 

•&(<_, o)— o 

Equations (2.19) and (2.20) with initial conditions (2.24) 
can be solved to obtain: 



f* _ X (T-TO 

^(X,T)= V, I e c Vf 

J o 

•&(*j x) = <f z f ‘e-*> (x ' x ‘' > yc£, 

<5 



(2.25) 



(2.26) 



It is customary to rewrite equation (2.26) in another form: 
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(2.27) 



-'o 

where 3£( T - -C ’) = 6 (2.28) 

The explicit form of the kernel (2.28) depends on the feed- 
back model under consideration [7] . Substituting equations 
(2.22), (2.25) and (2.27) into equation (2.18) leads to: 

-I- [(i- -i- CKCTT-'t'lv^C^/'T'^X* 

o 

+ v \ y i / e ' (2.29) 

J o 

which describes the increase in neutron flux above its 
steady value for one-group delayed neutron and Neufconian 
cooling model. 

The boundary and initial conditions for equation (2.29) 
are given by equations (2.23) and (2.24) 

° (2.30) 

q, ( *,o) = Fc*) (2.31) 

B. ASSUMPTIONS 

Before the analysis of the reactor may be undertaken, 
two further assumptions must be made concerning the feedback 
model . 

The first assumption states that the reactor temperature 
rises instantaneously with the flux (i.e., the prompt feed- 
back model). The kernel in equation (2.28) becomes 

[ 7 ] 
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(2.32) 



where 6(x) is the Dirac delta function, hence, the temper- 
ature difference in equation (2.27) can be written as 
follows : 

•'o "'o 

or T0'(*,'T) - (2.33) 

The second assumption states that the temperature 
coefficient of a reactor is negative. This is the case of 
negative feedback : 

|f < 0 (2.34) 

Therefore, the feedback coefficient a in equation (2.22) 

o 

is also negative: 

< 0 - ( 2 . 35 ) 

Inserting equation (2.32) into equation (2.29) provides: 





+ vV,/V V ’'' + 


(2.36) 


Let 


II 

1 

Q 


(2.37) 




0 

A 

« 

V 
1 

II 

O 

V 

J 


(2.38) 




vAAl = (3* 


(2.39) 
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Equation (2.36) yields: 




(2.40) 



and its initial-boundary conditions are given by equation 
(2.30) and (2.31). Equation (2.40) describes the increase 
of neutron flux above its steady value for a negative prompt 
feedback and one-group delayed neutron model, the solution 
of which will be established in the next section. 
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III. REDUCTION OF A NONLINEAR SYSTEM TO A LINFAR SYSTEM 



A. NONLINEAR TRANSFORMATION 

The Wilhelm's transformation [8] is considered: 






i - 



e 



(3. 1) 



or 






y (t . *c ) 



. )*(*, -t) 



(3.2) 



Since the neutron flux difference ip(x,r) is non-negative, it 
is logical to assume that: 

W'(*» < ’ £r .) '/f O -for" >1 i-r» 2) 3-net “~C O 

hence , 

i - ^ ( i - <z T ) y o (3.3) 



The nonlinear transformation (3.1) is generally valid, since 
it provides a unique interrelation between ip(x,T) >_ 0 and 
u(x,x) >_ 0, which is free from singularities for any point 
(x,t) . Differentiating equation (3.1) yields: 



_ 


1 Su. 1 __ &LL X ) 




(3.4) 


dz ~ 


\d% / 1- 




= 


K L ,Z 

e. V o 


2«(i-e K ‘)e K ‘( Vu ) 


(3.5) 
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Substituting equations (3.1), (3.4) and (3.5) into equation 

(2.40) provides: 



K. "C 

+ K (X U* & 


Q. V U, 


tt. **, 2 

, ^<(i-e )e (vu) 


^ 'r «r, 




Y *['- to - ^'^ u Y 



+ 



kt 

Ke u 



ztt 2_ 
oc C U 






Dividing both sides by 
give : 



+ V^ (r - V) e KT ^^V T l 



- e KT ')u(x,T') 



(3.6) 



KT 



™ y and simplifying 

[1 - £ (l-e KT )u] 2 



_ o* a^Ce -Q(W) 

K [‘“ k( , ' <2 ' CT ) u ] 



A* r, . - 0- e-) J Z fe' ^u^ } d rJ 



where 



Let 



and 



K 0 = K + > o 

K y O and 7l t > O 



V, = 



_ _ 1 V? u ) 



n _ A* r«-et(|- c* r ) ul f e I u 1 4 —' 

4 - r L ) J / I - ( , t «-) uU)T .) 



(3.7) 



(3.8) 



(3.9a) 



(3.9b) 



Equation (3.7) 

^ (i.-r) 

C>T 



becomes : 

= V 2 Ud<' r ’> + 7?, ( 5- ^ ) + T J 



(3.10) 
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The initial and boundary conditions are: 



^ £ e , x ) = 0 (3.11) 

= vp(i,o) _ F(< > ) (3.12) 



B. COMPARISON OF TERMS OF THE TRANSFORMED EQUATION 

Before solving equation (3.10) by successive approxima- 
tions method, it is necessary to impose some conditions on 
the terms and such that they are much smaller than 
the other terms in equation (3.10). A comparison of the 



terms of equation (3.6) indicates that the terms 



(W) 



K't 

6 and 



/ 2 \ KTJ KHP 

(V uje have the same order of magnitude as the term Kue 
Therefore, the nonlinear terms and are negligibly 
small in equation (3.10) for any time 0 < x < 00 if and only 



if the ratios 



Kue 



Km 



and 



Kue 



Km 



are small compared to 



unity, i.e., for x in D and 0 < x < 



JL 

« ue Kt 






« 



and 



(3.13) 



Vi 

Kua.*' r 




^Ci- 

te 



<r 



'u 



e 



ufx ,g)c<T* 






KT 






<< 1 



(3.14) 



Since u and Vu are regular, as Z tends to zero, the numer- 
ators of expressions (3.13) and (3.14) approach to zero, 
hence , 
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ku e' t ' t 



* o 



(■Lms 

r->o 






Kue. 



> ° 



For the time being, it is assumed that the conditions (3.13) 
and (3.14) are satisfied. They will be verified after the 
solution to equation (3.10) is obtained in the next section. 

C. LINEAR EQUATION TO BE SOLVED 

Let & = r], + % (3.15) 

Equation (3.10) becomes: 



* i ^7 ( * » 4 ) -+* £ 'T* } c 2. / ^3 


(3.16) 


J* 

II 

o 


(3.17) 


u(.*,o) _ F(X) 


(3.18) 



The transformed equation (3.16) now can be solved by the 
method of successive approximations. In this approach, the 
small nonlinear term '■&[t ,u (x , x) ] is treated as a perturba- 
tion. Hence, the initial-boundary-value problem is written 
as follows: 



nrrt 


w. + 




(3.19) 




0 




(3.20) 




F(i) 




(3.21) 
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where 






o 



and 



A = 1 }, + \= - fr x _* {i _ e * 



)u. 



(3.22) 



thus , the reduction of the nonlinear initial-boundary-value 
problem, equation (2.40) to a linear one, equation (3.19) is 
established by the nonlinear transformation, provided that 
the conditions (3.13) and (3.14) are satisfied. 
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IV. ANALYSIS OF A SLAB REACTOR CORE 



A. SOLUTION 

Consider an infinite homogeneous slab [Fig. 1] . The 
medium is of infinite extent in both y and z directions and 
of width 2a in the x direction (a is the dimensionless 
extrapolated half-thickness of the slab) ; the origin of the 
x axis is placed at the center of the slab. These specifi- 
cations reduce the spatial dependence of, the flux to 
variations in x alone. 




Therefore, equations (3.19), (3.20) and (3.21) become: 



T-'t ^x 2, * 

tk (*,<?) - F(x) 



(4.1) 
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where 



o 




and 





e U- U.-r'H-c 1 



(4.2) 

^ U— * / 3,. 



As seen in Section III, equation (4.1) will be solved by the 

l _ V_ 

method of successive approximations where the zero 1 " approx- 
imation corresponds to"£, = 0. 

1 . Zeroth Approximation (£, = 0 ) 

The zeroth approximation u Q (x,t) is defined by the 
following equations : 



^ LJ o 






dx* 


u 0 (± a,t) — 


o 


II 

y 


Fc*) 



Equation (4.3) is of the homogeneous, parabolic type that 
can be solved by either separation of variables or using 
the Green's function [9]. Throughout this work, the homo- 
geneous equations will be solved by separation of variables. 

Let U. 0 O,r) = Z(*)T(t) (4- 4 ) 



hence, ^ - XT' and ^ = X" T (4.5) 
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Substituting equations (4.5) into equation (4.3) and then 
dividing by XT yield: 

X' _ _X." = _ u Z (4.6) 

T ~ X P 

where y is the eigenvalue to be determined from the homo- 
geneous boundary conditions : 

S(ia) = o (4.7) 

Equations (4.6) can be written as 2 ordinary differential 
equations : 

X" + jU Z X = o (4.8) 

T' + y T = o (4.9) 

The solution to equation (4.8) is 

Z (x) = C, cos^rx -f (4.10) 

where c^ and are 2 arbitrary constants. The boundary 
conditions (4.7) give 

C ( ca.s^g +' = 0 (4.11) 

and C, - C 2 sin^a, = o (4.12) 

Subtracting equation (4.12) from (4.11) leads to 

=. o 

This implies = 0 or sinya = 0. 

If sinya = 0, then ya = n^, n = 0,1,2,... and equation 
(4.11) provides 

C, it = o 
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or 



C. = o 



and the general solution (4.10) is of the form: 



XC*) = C, sin 115 x 

^ a- 



(4.13) 



This solution is unacceptable since u Q (x,t.) is non-negative 
for - a < x < Si , hence , 



C * = o 



(4.14) 



and equation (4.11) gives 



C, — o 



Since c^ ^ 0 , then 



C oj p.a = 



o 



or the eigenvalues are defined by 



|X n = (2n + i)j£ ) rL = o,\,Z, 



(4.15) 



When equations (4.14) and (4.15) are inserted into equation 
(4.10) the solution of equation (4.8) is obtained as follows 



X (x) = C t cos (Zn + i)jL x 

Z a. 



(4.16) 



The solution to equation (4.9) is 



^ _(an + o*JL* 

T(r) = C^e 



(4.17) 



where c^ is an arbitrary constant. Combining equations 
(4.16) and (4.17) gives the general solution to equation (4,3) 
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°o 

= £ e» n e 



-Kn ^ 



cos^ n x 



(4.18) 



where are the eigenvalues given by expression (4.15). 

The coefficients B are determined from the initial condi- 

n 

tion: u q (x,o) = F (x) 



cx> 

o) = F(x) = ^ 3^0 os^. n - 



(4.19) 



r> = o 



The above equation indicates that the B n 's should be the 
coefficients of the Fourier's series expansion of function 
F (x) , i.e. , 



•*" ‘if 

-'-a 



(4.20) 



Replacing in equation (4.18) by its value in equation 
(4.20) yields 



u s Cx,~) 



=2 
n - o 



— r**' 

±e COCJ. c„x / F(f,)cos^ n ^di, 

J- a. 



- -y t: 

Since the series E e cosy x cosy £ is uniformly con- 

n=0 n n 

vergent for - a <_ x <_ a and 0 < x <_ <»[cf. appendix A] the 
interchange of signs E and / is possible, then, 



U Q (X.TT 



t} = f ^ 

Xa. L n=o J 



(4.21) 



Let 



s (*'&<•>•) = it e^hos^ 






(4.22) 
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Equation (4.21) can be written in the form: 



U o (X, "C ) 







(4.23) 



where G(x,£,t) is usually called the Green's function [8] 
associated with equation (4.3) 

2 . First approximation 

The first approximation is the solution to the 
following equation: 



> .f £ 

~ 1 
a.-r) - o 

w t (x,o) = Fc*) 



(4.24) 



where ^ is obtained from equation (4.2) by making i = 1 
and substituting equation (4.23) for u q (x,t). The solution 
of equation (4.24) is the superposition of two solutions, 
u and v to the following equations : 



£Ho _ (4.25) 

S Tt * 

u.(ta./0 = o 

U o(X,0) - F(x) 

u , = U a + V 

Equation (4.25) is identical 
the solution u q (x,t) is given by 



= 2!v + i ((X(T ) (4.26) 

v(fa.'c) a O 
o) = o 

( 4 - 27 ) 

to equation (4.3), hence 
expressions (4.18) or (4.23). 
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The solution to equation (4.26) is now assumed to be: 



=. ^ A n a)c o5^ n x (4.28) 

n -o 

the y n 's are determined by equation (4.15). The source term 
^(x,t) in equation (4.26) is then expanded in terms of the 
eigenfunctions, cosy n x: 

■£,(.*, ^) = &n iT) (4.2?) 

riTO 

where *£ n (T ) = JL f*t t ^ (4.30) 

J- a- 

Substitute equations (4.28) and (4.2?) into equation (4.26) 
to obtain: 



Y2 A n (T)COi^t n V = - Y2 pn A^JCc^X + '£Y' t)CC>& l Xn * 



n -o 



n- o 



dA 



where A (t) = . 

n at 



n 



or 



[ A n ( ‘ t) + f X n A n <T ) - t . V )] COS |X 



, x = 0 



(4.31) 



n =o 



Since the eigenfunctions, cos y^x are linearly independent, 
equality (4.31) implies that: 



A n ^> + ^ A n it) = -£/ CT) 



(4.32) 



n, — o / i , 2. , . • - • 



The initial condition is given by: 
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oe. 

v(x,o) = O = A n (o)co*^ n v 

n 5 o 

which is equivalent to: 

- O , rv = o,i/ -2-, (4.33) 

vS 

Multiplying equation (4.32) by e and integrating from 0 
to T yields the solution to this equation subject In the 
initial condition (4.33): 



A^c-r) _ e Hn f " d-t 

Jo 






where ^^(t) is defined by expression (4.30) 

- X 2U 



at 



or 



wi* 






(4.34) 



(4.35) 



rv = o, I, 2 .,. 3, 

Inserting equation (4.35) into expression (4.28) leads to 
the solution of equation (4.26) 






V^x,t') = 









n= o 



Since £.(£,t) is well-behaved for regular u(£,t), the uniform 

« -^(T-t) 

convergence of the series Z e cosy x cosy E , 

n=0 n 

[cf. appendix A] allows the interchange of the symbols of 
summation and integration in the above equation: 
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V(X,Tj - f f A (£'*') ^ ( r£ ^ 5 ^ 

Aj ^-a. ^ 



(4.36) 



'o •'-a 

From the definition of the Green's function, equation (4.22), 
equation (4.36) can be written as 



'O,'*) = ^ f GCx,$ ,~r- t) ek£<Lt 
Jo 



(4.37) 

— V. , ^ - -' U| v^ lV v - 

As seen from equation (4.27), the first approximation is 
given by 






f foe*,'} J T- t )&, (5, 

Jo -A a 



(4.38) 



. th 



. th 



3 . i Approximation 

It is noted that equation (4.1) for the i 1 " 11 approx- 
imation is similar to equation (4.24) for the first approx- 
imation with the same initial boundary conditions. 
Furthermore, the Green's function associated with a differ- 
ential equation only depends on the form of that equation 
and the boundary conditions [10]. Therefore, the solutions 
to all possible equations (4.1) with different (x,x) can 

be deduced from expression (4.38) by replacing (x,x) for 



£.]_ (x,t) 



( */ ~ o ^ ^ t ) •+ 



n e * ) J t l ci, t). 

i 



- 0,1,2,... 



(4.39) 



where (£ , t) is given by equations (4.2). Once u^(x,x) is 

evaluated from equation (4.39), the i^* 1 approximation of 
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the neutron flux difference will be calculated from the non- 



linear transformation (3.1): 






KT, , , 

e Ui ( x>to 



(4.40) 



»- £ (i-e^UsKii 



Thus, the analytical solution for an infinite slab reactor 
core with negative prompt feedback and one-group delayed 
neutron is found. 

At this point, it is necessary to verify the conditions 
(3.13) and (3.14) and the convergence of the sequence 
(u^(x,x)}. All these questions are discussed in the next 
subsections . 

B. VERIFICATION OF ASSUMPTIONS 

As shown in part III-B, the conditions (3.13) and (3.14) 

are justified at t = 0 and in its neighborhood provided that 

3 u. 

the function u(x,t) and its derivatives t — are continuous. 

d X 

These conditions will be verified for any x. 

For the sake of simplicity, the zeroth and the first 
approximations are used here to calculate the values of 
-£(x,t) and u(x,x) respectively. 

Expression (3.22) yields the value of (x,x) : 




where 




(4.41) 
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and 



Z f ^ -tTo ) 

7]' = A r,_ 2 < (.-e”)u d 1 e I tt(i^Ll±L ' 



(4.42) 



The fundamental mode of u (x,t) is given by equation (4.18) 

W e (x,T:) = & c e C> cojsu o x 



and so 



Let 



2id.‘0t-O = - & 0 ^e *5inp 4 K 

a/x-) = ^ e» b c<3 ^/ 



(4.43) 



The function a^ (x) is positively small for -a < x < a and of 

a -21 

order of =r B since a is very small (or order about 10 ) 



0 < &-, X zi 6 0 « > 

K 

Expression (4.41) and (4.42) are written as: 



(4.44) 



100' I Z*$>1 Mo 

IM- ^r. - -i*. v 



0c-2p*)t ~ i t i ‘ 



K^l-Q^oe" ^ A, c*>£ 



_ e J sin 

j 



(4.45) 



t = r- e h* **■*■ hh 

1 * 



(4.46) 



If the term a^e ° in equations (4.45) and (4.46) is neglec- 
ted for simplicity, it follows that: 
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I 111 « 



(4.47) 



2_ 3L C K "*^poJX 



><£i -f- a.cxJe* 



ck-hS )t 



J 



1^ 1 4 ['+«, 



4- A.CO-B 



CK-p, 



-K?m2 



f elZ!^' 



(4.48) 



a. Case K = U‘ 



Inequality (4.47) becomes: 



t -poT 



| rj‘ I / ^E>r KV / ^ e ° 

1 l( 1 *“ *c-«.3 K 1 



or 



r, 



Kue 



KT 



-KT 

/ £* 8 0 e 



2.«< £« 



-for ^ ^ o 



Assume a - 10 21 , B - 10 10 , K = 0.09, * = 10 -6 , then 

o o 



|_2L_ < 

| Ku 0** N 



-Zl \0 
2 < l o xio _. 

O . 0 4 



-9 



X \0 



or 



n, 



Kue 



XT 



^ ( -for ftnij 't > 0 



In the same manner, expression (4.48) is equivalent to 

- T 



3-, » 



* -K„T 2 (K„-pr;r' -*vrf ^ t -r' 

|r|^|4(J8 a a (i + Oji . ^ ^ (3 B. e 



o C 1 + A l ) 



$ Oie^-O < f- 

A ^ 



* n — A c r 



£. 



Ac 



e. e. 
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Replacing by (K+ ) in the above inequality yields 



i’ll i < *4- 



-KT 

e 



and 



'Au 



t Ic U ^ 



KT 



/ A* -KT *r —IC'T 

/ p C 6 = J£_ 

^ -X 









or 






KUC 



KT 



ft -for T o 



A c K 



Let 3*-10^, A ^ — 10 ^ and K - 0.09, then, 



tlL 



Kue 



l£T 



< T^f 



- 8 



( 

IT 



c . e 



KW & 



klT 



I < 1 ^ 



lo <0.0$ 

or cl>v<j T ; O < T < <> 3 



O .iO 



Hence, inequalities (4.49) and (4.50) lead to the following 
conclusion : 

I— 

I K U <2. 

This means that the nonlinear transformation (3.1) may be 
applicable in this case. 



« 



^or dnj OC. ■ T ^ oo 



C 4 • 5-' ) 
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b. Case K < y ‘ 



Inequality (4.47) can be rewritten as 



It I < ^ e 



2 



and so. 






KU Q, 






^4 T z 

</ So jJ-f ~v° r ^ ° 



Let y 2 = 0.09, K = 0.08, a = 10 -21 , B q - 10 10 , X c 10 6 
then, 



'n*. 



K UC 



K~t 



< 



- 1 / IO 

2x I O x IO x 0*0 9 



(o.oa) 1 



— 1 O 

3. x /O 



or 



K-U e 






4 ^ 1 -^ r an ^ T ^ ° 



(4.52) 



Similarly, inequality (4.48) turns out to be: 



- k/c (K.-uir)'t' *„ -K't r 1 

J ^ = " '] (4*53) 



Since, K - y 2 = K + X - y 2 = 0.08 + 10 6 
00 co 

it can be assumed that 

( K o - ) < 0 

then, inequality (4.53) yields: 



- 0.09 - -0.01 



n;i < i!&e w [,.^^ T ] < a" 

1 K m:-k. l j F. " 
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or 






XU e 






< 



e 



-(ZK- f^-o ) T 






2 ‘ Q 

If 2K - y_ 0, i.e., K >_ — — , it is obtained from the 



o — 

above inequality that 



_nk 



s< 



£ 



or 



Kuc^ » \ kty-l-Kc} 






|<U <2 



Kn: 



< 



for all T > 0 



i o 



- 8 



o. 08 [o. 61^ 



a 



x j o 



This implies that 



ni 



Kue 



Kt 



I «' 



for all T > 0 



(4 



Expressions (4.52) and (4.54) prove that the value of 
| i — j can be considered negligibly small compared to 



Kue 



y 0 2 

unity for any 0 <_ T <_ 00 if — K < y Q , 



c. Case K > y 

o 

It follows from inequality (4.47) that 



i ^ ( ic - 2 y* 

n: i < e 



for any T > 0 



hence , 



for anv T > 0 

'tcue^i N k 2 K - 



Let a = 10 21 , B = 10 10 , y 2 - 0.09, K « 0.093 

o o 



Then , 



.54) 
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JLL 



k uc 



or 



X 






< 

« 



-21 | o 

2x10 X I O XO-OC) 



- 8 ) 



( o. o 



— 2.1X10 



for any rp > 0 



(4.55) 



ku e 

To prove the ratio 



n. 



Kue 



K-n 



is very small compared to one for 



any 111 >_ 0 , it is divided in three different cases for the 
sake of simplicity: 



1 


T( (k-hZ)T 

If e 


« 


( i .C, 




(K- u» )t 

-e 


« 


1 < 


or 


m 


« 


1 /**/ K \ ^ _ 1 l„.( 0.0 9J>\ 




(ic-Hfl U3 0 ) d.oo3 W 2 !c/o'°) 




a 


« 


76(2.7 



Expression (4.48) yields the inequality: 



hi I < I s 



-Ko' 



J 






or 



and 



. O* n - Ko ^ CK,-p» a H 1 .*-(2 -Hot: 

’ll! < *-h e [ e - ‘J< tri^ 



| - 34*1 < - r for any t > 0 



Assume 6* = 10 8 , K = 0.093, K q - « 0.003 



36 



Then, 



Kue 



kt 



< 



-3 



10 



Z.- 



If a. e 



o. 0515 ^ 0 . 00 *, 
(K- Ho ^ 



-S’ 

A- * 10 << 



» 



for any t >_ 0 



c . e. 



1 



£ £7 



It is obtained from inequality (4.48) that: 



. I I * z *C«-Ho ^ ( 

< f *.e a,e 






S e 



< 7 ' r 






then. 



< ftp** 2. 

kA, 



1 ■n't | / £1» 

I (<u 1 He 2, A. k; a 



o ^=< 



K z A t 



K a A, 



-f or flny m ^ 



or 



I 01 1 I -8 10 -2.1 -i» 

1 £ / 1 0 * to x <o „ lo « I for any 

>Kue K ^\ (o.o^) z 0o" 6 ) 



3. _ 



If a, e 



Ck:- H©’)' 1 ' 



4. ■ £. 



K- 



v: 



\ ** (- ^ - 7U7 



1-4 ST 6 ) 



(4.57) 
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then 



I’ll | < 



— Kc 

e 



7 



&■ di: </ 



o'* rz " KoT C*.-H.7> T 

4 P ^o <2 






< 



4 (3.* 8 . 

it^T) 



It is deduced from the above inequality that: 




< 



-± .f* * Kaf < -*i * - 

X(K-C-^l) k (tc 0 -fjj) 



for 






^ o 



or 



n * j 

Kue </f I 



< 



-b 

lo 

o. o$3 y^.oa3 



-4 

l O << | 



(4.58) 



Therefore, inequalities (4.55), (4.56), (4.57) and (4.58) 

shows that the conditions (3.13) and (3.14) are satisfied 
for any x >_ 0 . In other words, for all cases, it can be 
written as follows : 



i 



kue^ 



« l 



for any x > 0 (4.59) 



C. CONVERGENCE OF THE SOLUTION 

Consider the sequence of successive approximations 
{vu (x,x) | u^ (x,x) is defined by equation (4.39)} (4.60) 



X r cU 



OkCx,T^ = u 0 O,t) -t J J G(x, %>, t- ( i, t) d-t 



(4.39) 



Whcrt- G> ( x t £, f x_ tr) — 



oC 



— ( -~ r ~ fc J 
2_ e COSpKCOi^J, 
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The Green's function G(x,£,x) is continuous in -a £ x,£ £ a 
for 0 < t <_ 00 [cf. appendix A] and bounded. Thus, the con- 
vergence of the sequence {u^(x,x)} reduces to that of 

{£• (x,x) } for 0 < t < 00 . 

1 . £i 

Since the ratio = 7 — is very small compared to unity 

Ku ± e KT 

for 0 <_ x £ 00 / we have: 

0 Kt 

-t. « k u. e 

For 0 <_ x <_ x^, the function u^(x,x) is analytically well- 
behaved, hence bounded by a number M and the inequality 
above can be written as 

K T , 

■i-l « KM <2 

which implies that the^^'s should be bounded by a quantity 

e ' 

-£• < %' -f« r ' o £ t: ^ rr , 

In addition, equation (4.39) gives the following recurrence 
formula : 

= f f (4 - 61) 

a - 

U ~ C> j I * 2. / * • * 

where the'^'s are evaluated from expression (4.2). Since 

— 21 — 8 

a is of order of magnitude about 10 and $* about 10 
[cf. appendix B] , the first term in the right hand side of 
equation (4.2) can be neglected compared to the second term, 

i .e . , 
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*• ■* J 0 . (i-e* T )u., ) fx,'t') 

It can be seen from equation (4.62) that: 

-ti >, ° 



(4.62) 



since the integrand under the integral is non-negative. 

The zeroth approximation yields = 0 and u q (x,t). 
Next, equation (4.61) gives: 



U, _ U 



■I 






-it 



where 

and 



G (x, § , <r_ t-J) >. o 

(.£,<=) >, o 



Therefore, we obtain 



- «o >, O 



or 



U 



(4.63) 



The values of corresponding to u q and u^ are given by 
equation (4.62) : 



t, * flTi . )u .] 2 f 

r L * J J o i-^0-e KT ')^„Lx,T') 



^ 0 U .] J 



^ _e 0 (T.V ) 



U, (x, tO^t' 



O £(ue KT, )u u.V) 



(4.64) 



(4.65) 



The values of u q and u^ do not much change the integrals 
in equations (4.64) and (4.65) but do change appreciably 
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the factors before the integrals. Furthermore, >_ u , 
hence , 

> t, (4.66) 

Combining equation (4.61) and inequality (4.66) leads to 
the following inequality: 

Uj- u, ^ o 



or u z ^.u., (4.67) 

Following the same procedure, it can be shown that: 



^0 


< 




< 


A 4 • 


• • • £ 




$ — 


(4.68) 


and u 0 




u, 


£ 


u i 4 • 


. . . < u- < 


^ t-V V 


^ . . . . 


(4.69) 



The sequence {*£^} is bounded and monotonic increasing, 
so it converges [11] . Therefore the convergence of the 
sequence {u^} is proved. 

D. ASYMPTOTIC EQUILIBRIUM STATES FOR AN INFINITE SLAB 
REACTOR 

Stability is one of fundamental problems in nuclear 
reactor design. Some questions are raised here. Is the 
system stable if the steady-state flux 4 > q (x) is perturbed 
locally or if the generalized buckling K is changed by 
changes in the absorption cross section due to control-rod 
movement? So, this sub-section is devoted to the investi- 
gation of all asymptotic equilibrium states (if any) from 
the solution obtained in paragraph A of this section. 
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* 









1 . Zeroth Approximation 



The zeroth approximation of the increase of neutron 
flux above its steady value is given by equations (4.18) 
and (4.40) : 






£ 



Ke 



-CKn-K)'' 



C :o5 









(4.70) 



r»s 0 



where the ' s are defined by equation (4.20). 

2 

a. Case K = = ( ^L) 

2 

All (y^ - K) are positive for n > 0 since 

Ko < < • * • <^n < Mn + , < • • • • 

Therefore, all exponential terms except the term in (y - K) 
will become negligibly small after sufficiently long time 
and equation (4.70) reduces to 



Sp X ^ <=)£ 't — S> (4.71) 



which describes the zeroth approximation for an asymptotic 
equilibrium state. 

For the case of no delayed neutron (8=0) and no feed- 
back (a = 0), expression (4.71) becomes: 

t°(x,*r) > B Cos JL x , 2Ls T — $> -4 (4.72 

This steady state is exactly the same as that obtained in 
Reactor Analysis [6]. 
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2 2 
b. Case y Q < K < y^ 

2 

All exponents (y - K) with n > 0 are positive, 

hence the corresponding terms in equation (4.100) decrease 

2 

with time t except the term in (y - K) which increases with 
2 

time since y Q - K is negative. Thus, as x goes to infinity, 
i^ o (x,t) approaches the following limit: 



^ 0 Cx,t) 



K ° ^ 









Since 



1c 






-C nf-Kl'c 
<2. COJ|U 0 X 



» 



cLS "T 






the stable equilibrium is 



> — 



(4.73) 



For the linear case, the reactor would be called supercriti- 
cal since the flux i^ o (x,x) will tend to infinity as seen in 
equation (4.73) by putting a = 0. Again, this result is 

consistent with previous Reactor Analysis [6] . 

2 

c. Case K< |X q 

All exponential terms in equation (4.70) decrease 
to zero as time x increases to infinity. Thus, the increase 
of neutron flux above its steady value is zero at the 
asymptotic equilibrium state for the nonlinear case. 

H'oC*' x ) > ° as t —?> <=& (4.74) 

This situation is called subcritical in a linear reactor. 
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2 . First Approximation 



The first approximation of the increase of neutron 
flux above its steady value is obtained by inserting 
equation (4.38) into equation (4.40) and putting i = 1: 



v h (Xf't ) = 



ACT . , , . ICT 

e ■+ A 



Lf7. 






• ) U 0 (x,t >- £ ( | _ <2 ")jf "J o, 0, I , T- O £ , f 4 , « <U- 



(4.75) 



where u q (x,t), G(x,£,x-t) and-£^(£,t) are defined by equa- 
tions (4.18), (4.22) and equation (4.2) respectively. On 

the other hand, -£.^(£,t) could be written in the form: 

IS.-O = + HiU.-t) (4.76) 



where o|(C/t) and are determined from equations 

(4.41) and (4.42). 



Only the fundamental mode is concerned with in the 

following calculations. Therefore, the values of u (£,t), 
3u o 

and G(x,£,x-t) are given by: 



Uofs.fr) 






f^ 6 (5,t) 


t 

= si^p 0 ^ 




-Mo(T-t) 


= -^e Co Y 



(4.77) 

(4.78) 

(4.79) 



Let I, 0 ,t) 



KT 



r -a. 



(4.80) 



tXiT) 




&(*>* , t-* )£, (s , t 



(4.81) 
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hence, the relationship between I^(x,t) and ^(x,!) is 
established , 

X, oe.-r ) = e* I z Cx.t> (4.82) 

and equation (4.75) becomes: 









q, UoC*, T ) -*- X, C^,-T) 

i - <sl C(_e K ' r )u o CK.T)-^l z cx f -n + ^ ; l t c^ ( -r) 



(4.83) 



Substituting equation (4.77) for u o (x,t) into equation (4.83) 
yields : 



MV*. - *) 



-ip. 1 - <•)'=' 

B 0 1 



(4.84) 



Taking the limit of the above equation, as x tends to 
infinity, gives: 



vp, t'i, t ) 

T -> •« 



Li rrv 
'T-? *« 



-(Mo-t-)T 

C OS}.*,, X 



x, 



(4.85) 






- & 



•K)T-| 

J cos V‘ 



x-si.1. 



six. 



It is apparent that the temporal behavior of i^^(x,x) depends 
on the generalized buckling K. The limits of 1^ and as 

X approaches infinity, are evaluated as follows: 



tlm. j f <3(x, £ ,T- Vjfc ( ( f ,f) <£b = 

^ T- 5 * 



r*L ( T 

-a J . 

(4.86) 



Thus, the problem is reduced to finding the value of the 
integral , 
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G( X, T.t ) ^ ( S , t- ) clt -far Hr t a r^g c. 



which is the sum of two integrals: 

f n ' r * 

J ^^/5 / T-fc)t| , (J | t3 cib and J 

6 O 

in view of equation (4.76). 

Equations (4.41), (4.77), (4.78) and equation (4.79) 

give : 



f ' , . _.,» t /»** -mH 



Let A, (v, 1, ) = “ 1^. 



(4.87) 



and ct,( § ) - 8c 



(4.88) 



The above equation can be written in the form: 



r - • -Ho 'tr' T -CH6-K)fe -Mo fc- 

eri’ olfc w A,d*.|te - *- dt (4.89) 



,(5)e 



Similarly, taking account of equations (4.42), (4.77), (4.78) 

and equation (4.79) yields: 



T -po^ 2 . (• 

J GT|^ cosM 0 * cas 



jri?* -> ‘U -‘^YfV^'^TUe 

e Q-»,e J j , „ 

-'o »-a,e r + A,e r * 
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Changing the order of integration in the above equation and 
letting 



& 0 cos f j -, y - cosZ P-3 



(4.90) 



leads to the following equation: 



f T I , a -<-{*S-Ko)x' r T cKo-K.^(r .u.'-f 2. 



a . Case K =JJ. q 

As t is sufficiently large, the integrand of the 
integral in the right hand side of equation (4.89) is equiv- 
alent to : 



-(p.,?-K)fc -u‘f 

_£ - £ 

a lG r^^a,e c ^'^ b “ 



» -*■ a,0&) 



(4.92) 



Substituting equation (4.92) into equation (4.89) gives: 



or 



J G r)J dt ~ A { Cx,l)e^° J 



cLt 



I-+ «•,(%) 



-UN 



G^J^t 2 i A 3 (x, S ) e X 



(4.93) 



where (y, 5 ) = _Ai C*> 1 ) 

\ °-i C ''s ) 

In the same manner, equation (4.91) can be written as 
follows : 



(4.94) 
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_ p 1 t r'*' -c^f- <»)-?' 

{ Gn 1 At- ~ A j_ C x , g )<2 . 0 + <a ( 

J '*• J 0 l-v^,(0 J_, 



2. , , 
,■) £. At At 



when t is sufficiently large or after integrating: 



f Grj 1 db ~ - A|OW)r T e_L 

-A ll Ho - L Ho - 



CH,r- ,<: o)'r 



K , 



+ 



H o - K 0 



] 



2. 

-Ho t: 
<2 



(4.95) 



But equation (3.8) is recalled: 





K a - K + X c 




or 


|<0- H^ = 




since 


K- = o 




hence , 




(4.96) 


Inserting 


equation (4.96) into equation 


(4.95) and putting 




=. Aa - (, ' vd '' 
^ C- 


(4.97) 



leads to the following result 

x PH 

or, as x is sufficiently large: 



G i} 1 At 



-H^T _A C T 

€ (-T - 



- A «' p N 

X- _ -i- \ 

X X . / 



r T; -Hot 

J <51']^ At — A + £ 'c 



(4.98) 



Now, adding equation (4.93) to equation (4.98) yields: 
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(4.99) 




C C -PX'* 

Stljic-t- Grj^t ~ (A 3 4-A^)€ x 

~ > €> J 6 



Integrating equation (4.99) with respect to £ from -a to a 
will give the integral I 2 (x,x) as defined by equation (4.81) 



,3. s'* 



I^tx.-c) = 



= J X A^Cxie 



-fN 



(4.100) 



where 



A 5 (>c) = ( {A i C*.5)-frA 4 (>‘/£0^i 

J _A 



(4.101) 



The integral I^(x,x) is evaluated by using relation (4.82) 
and equation (4.100) 






Iji. CX-zf ) 



KT -Mo ''f 

£ c A^- ( x) 't 



2 

Since K = y , then 
o 

I ( (x,x) t Ay (x) I (4.102) 

Replace equations (4.100) and (4.102) for I 2 (x,x) and 

2 

I^(x,x) in equation (4.85) with K = y o to obtain: 



Cirvv (*,-*) 



Li rrv 



l3» o C0f|Jc X + A^Cx) t 






Since the term A^(x)x is much larger than the other terms 
in the right hand side of the above equation for x large 
and A,, (x) is always positive [cf. appendix B] for -a < x < a, 
the limit of if>^(x,x) becomes: 



Lt'nv i^cx, x) 



Lll'V'V 

x ->•« 
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or I’.jyv 4> U,t) _ _!<_ (4.103) 

K 

i.e., iK (x,^) = — is the first approximation for a stable 

-L Ut 

equilibrium state when the generalized buckling K is equal 

2 

to the lowest eigenvalue, y Q , of equation (4.3). 

2 

b. Case K > y 

o 

2 2 
If (K - y ) is positive, so also (K - y ) is , 
o o o 

since K = K + X . 
o c 

For t large, the integrand under integral (4.89) 

becomes : 



e 

-CM 



“Ma" fc: 



-Cy.-vOt 

j <L r 






(4.104) 



Similarly, for the integrands in equation (4.91) 



-(Ho-<e)’ T ' 









'(lyf-K)T' 






Oc»- ic)-c' 






A iC?; 






(4.105) 



and 



-p^fc Cf**'*.)* _Z(M 0 -K^t 

e p-4, tl\e a-o^ido. J — Q. < V*) e 

Z O c - P* — At } t" 

— «-,U)€ 



(4.106) 



Substituting equation (4.104) back into equation (4.89) 
gives : 
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fo-V Ajc 


r 

~ A,(*,oe ^.t 


-wN 

- A fi (*, s ) e ' -X 


(4.107) 


J ‘ 

•'o 


J, 




where 


Ac (*<£ ) = 


A, (*, 5) 
«■, C 1 * ) 


(4.108) 


Inserting 


equations (4.105) and 


(4.106) into equation 


(4.91) 


yields the 


following equation: 








]® nL ^ - A * e A 

U C> 


— jM*)e 


cLtJiz' 



then, after integrating and some manipulations, it is 
obtained that: 












- X .1 

e i - e. 



*- A! 



1 



(4.109) 



The exponential terms in the bracket will approach zero as 
T tends to infinity. Therefore, equation (4.109) turns out 
to be : 



*'/s 






or 



r 



K-l 



(<- 



G ri 1 At ~ At. 2 . 

. 1 



Let 



A 7 (x,£) _ .ArC*^)^,fg ) 






(4.110) 



It follows that: 



J. 



r 

g Ab 



~ A ? (*. s )e 



(k - *f** )' t 



(4.111) 
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T 

Again, the value of integral / G£.,dt is obtained by summing 

o 

up equation (4.107) and equation (4.111): 



j Gt/t - ^ -f A 7 (*, 5 ) <? 



(K- Lp-t ) 



~ e 



CK- -(K-^)T 

I e -t + A 7 



(4.112) 



When t is sufficiently large, the first term in the bracket 

is negligibly small, compared to the second term since 
2 

(K - ii Q ) is positive, hence 



1. 



CK-ap* )~ 



- A 7 e 



(4.113) 



The value of ^(x,!) is deduced from integrating equation 
(4.113) with respect to £ from -a to a: 



(K- 2^o 



(4.114) 



Where 



X 4 - 1 6 

• ^6 (x1 = f A 7 (y ' fe *)^5 (4.115) 

5k 



A^(x,t) is defined by equations (4.110), (4.90) and equa- 

tion (4.88) . 

The value of I^(x,t) is determined from relation (4.82) 



X, ( X,T ) - 



V*' e 






(4.116) 



Inserting the values of ^(XjT) and I^(x,x) from equations 
(4.114) and (4.116) into equation (4.85) yields: 
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ilnru <1^ Cx, *t ) 



4 . ^pC 



Lc'rr^ 



{— 



B 0 e 






cx)e 



It 

C 0 SjU, 



+ A o «>oe 






* a,(«ie ' "- — a (x)e 



C/C- h o) r 






or 



tern- 
T — > <?o 



li rrt , 
T-> exa 



2(K "^ 1 )r r,7 , I 

_6r [^ 0 € Cosfi oX + A 8 J 






= Urn ( 



-(K-^^JT 

2> 0 c. CoSU 0 x -+ A g 



-u^yjx 
a,e +A,e 



-OC-pJ)t _CT 

' - ^ Age -f %Ag 



Since (K - y^) is positive, all exponential terms in 

the bracket go to zero as T approaches infinity. On the 

other hand, Ag (x) is always positive [cf. appendix B] for 

2 

-a < x < a, thus, the limit of ^(x,t) for K > y Q is: 



Lnrv = J<_ (4.117) 

T-$> *o 

which describes the asymptotic equilibrium state of the 
increase of neutron flux, 4>(x,t) - cf> o (x). This' is exactly 
the same result as obtained from equation (4.73) for the 
zeroth approximation. 

2 

c. Case K < y 

o 

The same procedure as carried out in parts (a) 
and (b) , the stable equilibrium state, will be found by 
using equations (4.85), (4.89) and equation (4.91). 






i - £ 



It t 



(4.118) 



" Kt 
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The above expression is obtained by dividing the numerator 

- (Vo-KJt 

and the denominator of the left hand side by e 
2 

Since y Q - K > 0, equation (4.118) becomes for t large: 



-cHr-Of -M-H 
•£ - £ 

I "Mo t -( 
i _ a.e. r " + A..e “ ' 






= -e 



-(p*- <)^ 



(4.119) 



hence, integral (4.119) can be evaluated as follows: 



sh 



jit — 



a, e | e <Lt 



A. e 



z 

-Mo T 



M*' K 



[>- e 3 



or 






M ^ (x, %ye 



-M 



2 _ 



(4.120) 



where 






r< 



(4.121) 



Furthermore, integral (4.91) is equivalent to 



fix 






<tt t A,£ 



A -fit r CHa-Koyr 
■ — 1 A ^ g. j 2. — I _ T 

Mo' 2 '- ^ Ho 1 - Ko 



It is assumed that: 



fV - > 



then, for t large: 
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-K(fr 

A z e. r 

(nt - K.) 2 



' Ko ) ^ 

-e 



-Ko'C 



— A lo C*, ^ ) e 



(4.122) 



where 



A, 0 C*.-?) = A »^«) 



(4.123) 



Adding equation (4.120) to equation (4.122) yields 

2 _ 

“ Aq ( x ,'%)(2 4 - A tc C K ,%)£ 



& 



Since > K, 



thus , 



J" Gi, si A, 0 (x. 



- Kot 



(4.124) 



and 



3 2 (v,-t) = [ c; A m cx)^.' 

“ 4 . 'J O 



_k.T 



(4.125) 



where 



A„ Cx ) 




t«.S) ^5 



'■'-A 



(4.126) 



and 






K-C 

e LchaI ss 



Oe-tXo)T 



— A M o;e. 



a,t 



(4.127) 



Finally, inserting the values of 1^ (x,x) and I-^(x,x) from 
equation (*.125) and (4.127) into equation (4.85) yields the 
first approximation of the asymptotic equilibrium state: 
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Um, if/ i*,-?) 

T^> oo 




B. 



I - A,c*)£ °-v- <x.c + )a 



-tp^-K.)T _^ C T 





•Ot —k 0 t —X,x 

- •£ A „°‘>e + iA „«»< ‘ 



Since all exponential terms decrease as t increases, then: 

Li »rw tp’i (x.,t) = o (4.128) 

T— 3> PO 

Thus, the same asymptotic state as for the zeroth approxima- 
tion is obtained. The same conclusion can be made even if 

2 2 

K q is assumed to be greater than y , instead of K q < y o « 

E. PHYSICAL INTERPRETATION OF RESULTS 

It can be seen from the zeroth and first approximation 
that higher approximations are not required for asymptotic 
equilibrium states in a nonlinear system with negative 
prompt feedback and one-group delayed neutron. The stable 
equilibrium states in the first approximation are verified 
by the results of Kastenberg and Chambre's stability 
analysis [1] . Furthermore, in the absence of feedback, the 
results obtained here reduce to those of the linear reactor 
analysis as mentioned in part D-l. 

The equilibrium states are established in a nonlinear 
nuclear reactor at different levels which depend upon the 
value of the generalized buckling K. The results are 
interpreted as follows. 

For t < 0, the reactor was in steady-state with flux 
<P Q (x) . 
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At x = 0 , the flux is raised to an amount of F (x) from 
an external source. What happens to the flux distribution 
after some time (t > 0)? 

TT 2 

1. If K > / the neutron flux will reach the stable 

equilibrium state of value <£(x) = <j> (x) + — after the 

transient given by equations (4.39) and (4.40) has subsided. 

TT 2 

2. If K < (■=— ) , the neutron flux will return to the 

Z 3 . 

initial steady value <f> o (x) as equilibrium state. 

On the other hand, the equilibrium state is weakly 
dependent of delayed neutrons since the fraction of delayed 
neutrons (3) is very small compared to unity in the expres- 
sion of parameter a: 

©< = lO-f&Wa** | 

If delayed neutrons are neglected (3* = 0), equation 
(2.40) becomes 

^ 'iL ( *, t ) — (x,t) -t '■Vc k,t ) - <P*7x .t - ) (4.129) 

d T ^ 

7T 2 

For K > (~—) , theorme I of Kastenberg and Chambre [1] gives 

the asymptotic stable equilibrium state of the system 
(4.129), which is the positive (nontrivial) solution of 
equation : 

"Ciryv [" K<4>(*.'-r'>--c*iy 2 7x / -c)1 _ q (4,130) 

o L 3 x 1 > 

Y 

It is apparent that \p (x) = — satisfies this equation, thus 

Gt 

is a positive solution to equation (4.130). Therefore, the 
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result obtained in part 1) is verified by this theorem for 
it 2 

K > i^) • Similarly, theorem II of Kastenberg and Chambre 

77 2 

shows that the equilibrium state for K < (■=— ) is \p (x) = 0; 

i.e., the neutron flux returns to the initial steady state 

<P o (x). This is again consistent with the result in part 2). 

2 

Theorem III states that if K = y o , the asymptotical 

stable equilibrium of the system above is also the zero 

state (ip (x) = 0). This completely differs from the present 

result as described in part 1) . From the physical point of 

view, if the reactor eigenvalue is not changed, then the 

perturbed flux should eventually return to the initial 

2 

equilibrium state \p (x) = 0. The fact that for K = y Q , 

the first approximation contains a secular term at large 

r [cf. equation (4.102] indicates that this is a special 

case for which the Wilhelm's transformation may not be valid. 

To draw a definite conclusion at this point, it appears that 

2 

the first approximation of the solution for the. case K = y 
requires further considerations. 



58 



V. ANALYSIS OF A FINITE CYLINDRICAL REACTOR CORE 



A. SOLUTION 

Another geometric configuration which is of practical 
importance in Reactor Design is the cylinder of finite 
length. Consider therefore a bare cylinder of dimensionless 
extrapolated radius R and height 2%, as shown in Fig. 5.1, 
with the 2 direction along the axis of the cylinder and the 
origin of the coordinate system located at the center. The 
flux in this reactor depends upon r,z (cylindrical coordi- 
nates) and t. The reactor equation is now [cf. equation 
(2.40)] 



2£(r,*.r) = + K 



cx ' 2 -*- (S* 




(5.1) 



o 



for 



0 4 r <£ K 



where 



VV - ' A. 2. frl±) -v Elf 
r ^ r \ ->r I 



(5.2) 







Figure 5.1. A Bare Cylindrical Reactor. 
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The initial and boundary conditions are: 
vf-Cr, e , o) - F tr,*) 



Lr -finite. 



(5.3) 



Cr, ± l,t) = ° 

Using the same transformation (3.1), equations (5.1) and 
(5.3) are transformed to linear equations as follows 
[cf. equation (3.31)]: 



r) _ ^ u.j cr, e , t ) ■+ ^ C r, z , -c ) (5.4) 

u iCt, e, o ) = FCr,^ ) 

® 

(5.5) 

U ; (0, e, T ) -finite. 

u l (r, ± / , T ) =• <5 

where -£<> = ° (5.6) 

and 






f ['«■ 



■ 0-e 



f "t _Ko(T-t'| . , 

:r. u 1 Z e U;_|(t*< T ) <*T (5.7) 



The methods used in the analysis of the infinite slab are 
easily extended to the analysis of the finite-cylinder 
reactor. The procedure is straight forward and involves the 
use of the technique of successive approximations and separa- 
tion of variables. 
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1. Zeroth Approximation (£ = 0) 

The approDriate form of the transformed equation 
which gives the zeroth approximation is: 



— <=“ _ JL £_ (r—'\ 
vn: r \ a r / 






3 e ■ 



(5.8) 



u oCr, -t, o ) = P (r,*) 

= O 

U 0 ( o, i. T ) is iCrt'vte- 

u o (r , ± l , T ) - O 



(5.9) 



A solution of this equation is assumed to be in the form: 

U 5 (r, a, t ) = •X(r') Z(e) / T(r) (5.10) 

The substitution of this expression into equation (5.8) 
yields : 



T <t\ 



T(t) rjjto 



4- (rT) = j-Lll = - M- 3 
drK ' Z (*) 



(5.11) 



where U is a constant to be determined from the boundary 
condition (5.9) or 



£(-0 = Zu) = o 

It follows from equation (5.11) that: 

Z% n*Z = O 



(5.12) 



(5.13) 
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and 



— -t f> 1 ' = — 4- 

T ‘ rX 4r 



(rX'J = - 't 



(5.14) 



is another eigenvalue to be determined from the boundary 
conditions . 



X ( e . ) = o 

X ( o ) is -f irvlte 

Equation (5.14) can be written in 2 different equations 



(5.15) 



and 


h ( rX ' ) + ^ r:E 


= 0 


(5.16) 




T ' + (jaV'Ot 


= 0 


(5.17) 


Solving 


equations (5.12) and 


(5.13) gives: 






p-rv 


= (2rn-i)22- . 'tt-Oj 

24 J 


(5.18) 


and 


£(*) 


= . C. 


(5.19) 



Equation (5.16) is a Bessel's equation, the solution of 
which is of the form: 



X(r) = Cj 0 (WvJ) + PY 6 C'tr) 
where C,D are 2 constants 



( 5 . 20 ) 



5C(o) = cj 0 to) 4 - s> Vo ( o) 
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3M! 

H v 



hence, D = 0 since y (0) = 00 and X(0) should he finite 
[cf. equation (5.15)]. Furthermore, the boundary condition 
X(R) = 0, yields: 

J e (-Tr) = O (5.21) 

If oT 2 ' “ ' * ‘ are P os i t; *- Ve r °ots of equation 

(5.21), the solutions to equation (5.16) are given by equa- 
tion (5.20) 



XCr) - C J 0 (tn r) 



and the solutions to equation (5.17) are 



TO) 



E e 






(5.22) 



(5.23) 



where E is any constant. 

Combining the product of solutions (5.19), (5.22) and 

(5.23) leads to the solution of equation (5.8) 

= £ A 

m,n r o 

The constants A are defined by the initial condition 

m , n 

=F Lr,*) = £[£ 

n:« J 



e 



i -*■ ^ 






(5.24) 



This equation implies that: 



A. 



i 

..nXC'C*'*) = x FCr*,f )Cosp n jc = fv,Cr; 

J-£ 



( S'. 2.0 



rv =. O , 1,2., 
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Expression (5.25) is known as the Fourier-Bessel series 

representation of F n (r). Since are the positive roots of 

equation (5.21), the coefficients A should have the 

m,n 

form [12] 



i, n = 















or 



A 



rh , r\ = 



J o J-i 



f (5 ' 26) 

J 0 J-i 



Substitution of expression (5.26) into equation (5.24) gives 

the solution u (r,z,x): 
o 



U 0 fo*/T) = 



fa • 5 ) fas-d < 5 • 2 7 > 

L J.VnA) J ( 



y o - i. 



Let 



G(<\* ; »1 , $,t) = 2J1 J’ 0 (X t ()ccs^<u>s.^l (5.28) 

u ^ 



rn , r» = a 



G (r , z ,n , ? , t ) is called the Green's function of equation (5.8). 
Then, equation (5.27) becomes 





f * f 1 


U oC*\ T ) 


J 6 J-i. 



2 . First Approximation 

The first approximation satisfies the following 
equation : 
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7 t 



_L 2. fr 1 
- r 3r^ ar j 



4 



9 



£, C r / * 



(5.30) 



(5.31) 



C r, i, o) - F (r e) 
“,U , 2 ,r) = O 

U , ( o, a , T ) LS f in4e 

U, «,T) = o 



where ^ (r,z,T) is determined from equation (5.7) by making 
i = 1. 

The solution of equation (5.20) is the sum of the solu- 
tion u Q to equation (5.8) and the solution v to equations: 



2 y _ J_ 2 .fr?X.) 

a-r r ^rv. / ?r/ 


. i!y + f 


( 5 . 32 ) 


< 

P 

o 

II 


o 






O 




v (o, Z., -T ) Li -Rin 


tie 


( 5 . 33 ) 


vC r /i£,T) = 


o 




u = 

1 


u. + v 


( 5 . 34 ) 


The general solution 


to equatoin ( 5 . 32 ) 


is of the form: 








vm,' o = 2^ J A p ) co - s i u ^ 2 


( 5 . 35 ) 



r» ,n = o 
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Expanding-^ in terms of J q (y^r) cosy^z yields: 
( r ' * 1 T ) = 'll C m,K f ' t ' , Js('^m r )^^ A 



(5.36) 



M s o rv» s: o 



where E C J r) are the Fourier's coefficients of the 
m=0 m ’ n ° m 

above series, i.e.. 



( 

Z- = -1 J £, C r, j, , T ) | 



This equation implies that C (~) are the Bessel's series 

m,n 

coefficients : 



C*..rv< T > = 



.e ,< 



(‘Vn» R ) -x 






(5.37) 



Taking derivatives of equation (5.35) gives 



'tL 

9v 

'"d r 



3 2 v 



9£v 

9^ 



CHJ 

\ C) Co XU S' 

z_/ a.t * " 

m |fl? o 

- X2 ^ ^ c ° s \ x ~* r 

D 

- £ 6„,.p.(V)- ±- • J, C-^r) ] 

fv * 

■*» » Vi * 0 
■o 

~ ^ Pm &>!"•,« C-oSy^-Z: 



(5.38) 



m, * » o 



Insert equations (5.38) and (5.36) into equation (5.32) to 
obtain: 
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o 



CXJ 

L [#-" C m .„ 

m,nsO L 



This equation is true for any 0 £ r < R and -% < z < l if 
and only if: 



4 *=- -v 

4 T 



rrc. *o = 






(5.39) 



The initial conditions are given by: 



oO 

v(r,*,o) = O = ^ C*) J 0 (^r)cos^ n A- 



Dr\i h, e ® 



or 



£> (o) = o 

ID 1 1\ ' 



(5.40) 



The solution to equation (5.39) can be written as: 






Jut' 



where D are constants determined from the initial condi- 
m , n 

tions (5.40) 






Therefore, it is deduced from the above equations that 



Bvn.n ( -' c) 



J 



<*■*= 



Replacing equation (5.37) for C (t) leads to 

m f n 
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t> It) 

rh , n 



Te5WT) J j ~ (5 - 41) 



Then, it follows from equations (5.35) and (5.41) that 



r R r i 



v(r, js-,t) = 



•"X m,». ■=■ o 






^ J 4 CC -l) M „* COJT p „ 5 fe <tn <l<r 



3, 



or, in terms of the Green's function (5.28), the above 
expression becomes: 







f. r t - 

'£,(' t ), 5 , 1 ) <=>(>, /^-0 

* 



(5.42) 



It is noted that this form of solution is analogous to that 
one derived in equation (4.37) for a slab reactor. 

Now, combining equations (5.29) and (5.42) yields the 
first approximation: 




(5.43) 



3 . ith Approximation 

Since equation (5.4) is similar to equation (5.30) 
with the same initial boundary conditions, the ith approx- 
imation which is the solution to equation (5.4), can be 
obtained from expression (5.43) by substituting (n , £ , t) 
for : 
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(5.44) 




Where ^ (n , £ , t) is defined by equation (5.7). Equations 
(5.44) and (5.7) will give all possible successive approxima- 
tions of the solution to equation (5.4) provided that = 0. 

Finally, the ith approximation of the increase of neutron 
flux above its steady value is calculated from the transfor- 
mation (3.1) 

ICT’ 

=r ^ — '±UL±1^L ) (5.45) 

1 - £ (l- e K:T ) r,3r ,-c ) 

Thus, the analytical solution for a finite-cylinder reactor 
core with negative prompt feedback and one-group delayed 
neutron is established. The convergence of the series 
{u^(r,z,x)} can be shown by the same method as indicated 
in part IV-C for a slab reactor. 

B. DISCUSSION OF RESULTS 

From the solution (5.44) and the transformation (5.45), 
the asymptotic stable equilibrium states can be found at 
different levels by the same approach as applied in a slab 
reactor. The proof is not repeated here but only results 
are summarized as follows. 

a. If K - v 2 a ttj 

After a long time, the increase in neutron flux 
above its steady value will be constant with time and reach 
the value : • 
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( r, t , x ) 



(5.46) 



H2 cc,J r° a 

l+ ^ A o 7 4 ( Y . r )C“^.^ 

where A q is given by equation (5.26) 




2 . 




ri 

r)X(^r]jF(q / | 



b. If K > y 2 + gr 2 

o o 

The zeroth approximation for asymptotic equili- 
brium state is : 

C r / ■* » 00 ^ = — (5.47) 

c. If K < y 2 +r 2 

o o 






(5.48) 



It means that after the transient the neutron flux will 
return to the initial steady value as the stable equilibrium 
state. It is noted that in the absence of delayed neutron 
(8 = 0) and feedback (a = 0), equation (5.46) becomes: 

H'oCr, *,*) = A o^(X r ) coS ^o Sr (5.49) 

for K = y 2 + V 2 . 

o o 

This result is exactly the same as that one obtained in 
previous work [6] . 

2 . First Approximation 
a. If K > y 2 + 

The neutron flux will approach the stable equili- 
brium value: 
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(5.50) 



0(r, */ »° ) = t 0 (r, 2 ) JL 

as time goes to infinity. ' 

b. If K < y 2 + IT 2 
o o 

The equilibrium state of neutron flux is the 
initial steady value 4> Q (r,z) 

^(r »*/ r , (5.51) 

Therefore, the same physical interpretation can be made 
for a finite cylinder-reactor core as for an infinite 
slab-reactor core. It seems that the asymptotic states for 
a nonlinear nuclear reactor do not depend upon geometric 
configurations. This is true since the transformation is 
independent of coordinate systems and the transformed 
equation is maintained in the same form. 
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VI. CONCLUSIONS AND RFCOMMENDATIONS 



The objective of this thesis has been to provide an 
approach for directly solving the nonlinear space-time 
reactor kinetics equations. The theory is developed for an 
one-velocity, bare, homogeneous, one-group delayed neutron 
and Neutonian cooling system where the temperature rises 
instantaneously with the neutron flux and the temperature 
coefficient is negative (negative prompt feedback model) . 

The fundamental tool in this analysis is the Wilhelm's 
nonlinear transformation [8] that transforms a nonlinear 
equation to a weakly nonlinear one in which the nonlinear 
term represents a small perturbation. Then, the reduction 
of the later equation to a linear, parabolic type is estab- 
lished by any method of approximation. The successive 
approximations technique has been used in this thesis since 
it yields a general form of the ith approximation for the 
solution and it is easy to simulate the problem by digital 
computers if need be. 

Furthermore, this transformation is independent of the 
system of coordinates. In other words, the transformed 
equation is maintained in a general form whatever coordinates 
are chosen. Therefore, the same conclusions of stability 
can be drawn from the resulting solutions for reactors of 
different geometric configurations. The asymptotic stable 
equilibrium states are governed by the relationship of the 
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2 2 2 

generalized buckling K to m q (or y Q + 'f ' ) , the lowest eigen- 
value of the associated linear Helmholtz's equation. The 

K 

value, <j> (x) + — , is the stable state of neutron flux when 

2 2 2 2 

K > y (or y + T ) . When K < y , the neutron flux returns 
— o o o o 

to its initial steady value <j> (x) after some transient. 

Then, these results are compared with those derived from 

Theorems I, II and III of Kastenberg and Chambre (1] in the 

case of no-delayed neutron. As mentioned in part IV-E, the 

only inconsistency of the present work with Kastenberg and 

2 

Chambre' s is the stable equilibrium state for K = y Q . But 

if the comparison is made with linear reactor analysis (no 

delayed neutron and no feedback) , the zeroth approximation 

(cf. see part IV-D-1) will give an asymptotic stable value 

which matches very well with the previous result [6] . So, 

for future work, it is recommended that a verification of 

2 

the first approximation at K = y Q by another approach, say 
digital computers , be made to insure its consistency with 
the theorem III of Kastenberg and with the physical point 
of view. 

It is noted that the successive approximations method 
can be applied directly to equation (2.40) to obtain the 
solution provided that the coefficients a and 6* should be 
very small compared to other terms. On the other hand, the 
Wilhelm's transformation reduces equation (2.40) to equation 
(3.10) where the terms n 2 and n 2 become smaller, even if a 
and 3* are not appreciably small. 

Finally, it is assumed from the beginning of this work, 
that the feedback coefficient a Q is negative. May the 
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transformation work with a positive feedback system (a > 0 

or a < 0)? Moreover, if ip >, 0 r u(x,t) should be non-negative 

K 



and less than the value, - 



a(e -1) 



Otherwise, the func- 



tion ^(x,x) will be negative, i.e., the range of u(x,r) is 
limited in this case 



o $ ct < _ 



K 



( e * x -0 



o' < O 



and the singularity may occur when u(x,r) is equal to 
K 



a(e KT -l) 



at any point (x,r). Therefore, it is not con- 



venient to apply this transformation to a problem with 
positive feedback. 
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APPENDIX A 



« -y 't 
n 



UNIFORM CONVERGENCE OF THE SERIES E e cosy x cosy £ 

n n n 

n=0 



TT 

Consider the above series where y = (2n + 1) ~ — 

n 2a 

will prove that this series converges uniformly for 
-a <_ x, £ <_ a and 0 < x q <_ x £ <». 



We 



✓ \ i ~y» r i , -r* '* / -r*> ■■ , 

(») |e | <; e. e x fir 



. -pi ' r « 



- 



Tp 4 C ao 



y*x Q 

Expanding e into Taylor's series yields: 






= I- 



-V JL (>* T.)^ 
• 



hence , 



M^To v 1 

e " > KT. 



and 



-For, 

e .4 



4A 



pi T ° 



7T 2, To^nti)' 



Therefore, inequality (1) can be written as 



, -Un? I ^ a, . /-«<*., t 

(2-) | e coi^xcoc^? 1 < ~~i~x 7— T Z. ~^ r | 

TI / 0 + 1; £ c 



- a < *, j: C ^ 

cx 3 



It is known that the series E 



is convergent, then 



n=0 (2n+l) 



2 

* -y t 

inequality (2) implies that the series E e cosy x cosy £ 

O n n 

should converge uniformly for -a <_ x, £ <_ a and x q <_ x <_ 00 . 
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Let 



» x cc,s b" ? 



£ C * • s , ^ e <^ s f* 

n = O 

Since each term of the series above is continuous, the func- 
tion S(x,£/X) is also continuous and bounded [11] for 
-a £ x, 5 f. a and t q <_ t <_ «. 
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APPENDIX B 



1. Calculation of the Coefficient (x) in Equation (4.100) 



Assume K 






(— ) 2 
' 2a 



From equation (4.101); 



A 5- (V ) = J Aj (*, £.) dj* 4- ^ A ^ (>c, 5 ) ^ 



where 



^ (*'5) 



A , C »> ?) 



2<x 

«-i< 




c.osHe*c o SFoi;Sin 



I -t 






( 1 ) 



and 



AiO±£i> = _£? K o»su 6 xras , 'p 6 S^n-^-g. 0 ^L { . 6 |] 
n K 



Let us evaluate integral / A^fXfOd? 



-a 



f. 

f = - — s' f *- 1 jj 



p CL ^ 

< *‘ C ^ 6 1 -t ^ 



Let 



P 

3 



CoS |^ 0 £ 



L <^ = - sin p t £ 



( 2 ) 



(3) 

(4) 



Equation (3) becomes: 
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( 5 ) 



The integral in the right hand side of equation (5) is 
separated into 2 different integrals as follows: 



But the first integral in the right hand side of equation 
(6) can be calculated from C.R.C. tables, 

j = x (7) 

The second integral can be evaluated after changing of 
variable: 

Let ■ X. - vj -+• p 






hence 




(8) 



Since it is obtained from any integral table that: 
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f )Ll £ ±2f^±2l = + b f f(x. 

J x J V ^V^.-p* 



:+(/-P^J 






^^-X. l +i]?Xtl- p 1 " 



( 9 ) 



and 



f - ^ = sln-'cx-,,) 

J l/— p" 2 - ^ V^f ^ 



sirT’cj. 



( 10 ) 



and 



l 



^ _ _l fin - 1 ( P^-*-' ~ 

-X V^Z%UfJC*-\- p2- \^7 \ a 



sin"' fei±±A 
V Vf + P/ 



(11) 



\/ pK-i 

Substitute equation (10) and equation (11) into equation (9) 
to obtain: 



\Z-xV^>cf i- y>*- _ Vi- -pp5>rT’^_ x/p u -i (12) 



which is integral (8) . Thus we have 



f It— lu, - - 1 +■ j-l/ r p u -i + \/F^r c> n - ' -j- 

*■'6 ‘3-+ b ^ r 



K 



Since p = is very large, hence, 



(13) 



^ « 1 



and Vp*.i = t( 1 1 - P 



J 






5 m~ ^ _L ^ 

P P 6 p 2. 



Equation (13) can be written as follows: 
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1 



TL. _ _J_ 

4p Sp 2 ' 



(14) 



Now, inserting equation (7) and equation (14) into equation 
(6) yields: 



I 






& 0 



£1 . JLu = l£.-bfiE.__L_ > \=_i- = 
c y- + p ^ 3p*-/ 2-p 3K 



which is the value of the integral in the right hand side of 
equation (5). Thus, 



J A 3<>, O A 5 = " 

- 



4- i. ^ 

3<jt< COS Po k 



But K 0 = 2a 



Then 



J 



(15) 



Next, evaluate integral / A. (x,£)d£ where A.(x,£) is given 

-a 

by equation (2) . 



J' 



A q(*, -g) d? = 8 a c°Sfj„x. f( |_+ ^ 

^ C <J_ 



= fr f u; 
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Let y - p 0 c - 53 - 5 



dy = ^ d 5 



hence , 



f ^ = z r co - s ^6 $ ^ 



- ** 
~TC 



r ?: 



Cos 



‘a “V 






__ if A. 



( jr \ ~ 

~IL [ *t ' 



<SU 



(17) 



and similarly, 




8 a (18) 

3 TO 



Substituting equations (17) and (18) into equation (16) gives 




A 



CCC I 






t> + f 






1C 






* r -f* + 



L- 



6 



3 7Ck V 



L> ] 






Gfi S 



t 1 



(19) 



Thus, Aj.(x) is obtained by summing up equations (15) and (19) 





~ t &0 CoSjJ*) c + 



i- Lz \ £> 0 x 

VAt, 3 TK. A 0 / 
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or 




L 

3 TTKAc 



k_ X 
3 a- 1 - 



S„ COS^ e X 



( 20 ) 



Let us consider the order of magnitude of each term in brac- 
kets. To do this, we assume the following data [13] 

a £7 S> - O.oo 3 3 > 



^ =■ O. * 

^ = d>. 0(3^7 $ / cm'* 

-y 

V z. H-- I x to <-r* fitc. 

A = <?• 8 . 2 .^- 




0 . 0 0 £“ 2. CkV! 



dj, = 3- 07 ? c^- / -g 19 ^ 

-T Q .1 

* e ^ - rx/o 0 ic 
-e = 2.00 = 2aoK3.S2.mo' i+ '£ 

%-is.i i <>•>» j^rssion 



Therefore, we obtain: 



K = ( 


: s-j 1 = 1 















K -fe 0 -i ~6 0 = I. 



ll^ - 4 ^ £ - U k 10 cm 1 " ate *< 

^ — | i i-^.) c *i>~ Y -j. | a I =: S“. ■T'x./o cm 1 o*. o 



*c. = 



o. 8 2-<+ 






&£l%>o 7 ^0 -oori- 



— 3 . 7 . x 1 0 



— < 



Cci 



t m t>v-» j|’on 






L sL 

i *TC K 



-2o 

4*- “7 >•/ O c 



*- 

rvw/ /Occ^ . 
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O * o O 3 6 * 2 — 



( 21 ) 



4 ^ — - (~> fa ~ $ fe-c 



3 



jfk _ / 6. \/£M = 

TKA7 " (-3 TT^JUc ) 



- 2 *- «L 

/ « ~7 */ O Cm .s «- c- 



*2- "Tf p( ; . -» ^ 'V. 

— — . — =: M- - & xlD cm o^c^ 

3 K 

^ _ Q gj 

^ — o ■ oe> %C"i- A c- = o -o & 'iG'i-'/. 1 -2. *■ l & — / . f x / o 

Thus, comparing equation (21) to equations (22) and 
yields : 



£T 

A c 

€. 

Ac 

and the value of A^ (x) is approximately equal to 

Av Lx) — Jrr- S’, £&£ 2L. X 

Hence, A,, (x) is finite and positive for -a < x < a. 

2. Calculation of the Coefficient Ag (x) in Equation 

2 

Assume K > y 

o 

From equation (4.115): 



1 

3 



7TK A(_ 






and 



— JL. sL 
3 <t v k 



<C 



A 



a <*> 




( 22 ) 



(23) 



(23) 



(24) 



(4.114) 
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-A z ±i 



coc^i fto 5 ) 



where 




C<~ h-o 1 ') 



$2* 

°-<^c. LK-^o 2 -) 



y /"» V 

A<j(>0 = -A * g ° or,u 6 x /oc 

^IcAcCk-^o 1 -; 1 J _* 1 { 

Since the value of the integral in the right hand side of 
expression (25) is given by equation (18), hence. 



A 8 M 



&_ 

3 



- -&* ** g J oac X x 
KAcCfc- po 1 ") 2 - <a 



(26) 



The equation above shows that Ag (x) is finite and positive 
for -a < x < a. 
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